On the Causality and Stability of the Relativistic Diffusion Equation 
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This paper examines the mathematical properties of the 
relativistic diffusion equation. The peculiar solution which 
Hiscock and Lindblom identified as an instability is shown 
to emerge from an ill-posed initial value problem. These do 
not meet the mathematical conditions required for realistic 
physical problems and can not serve as an argument against 
the relativistic hydrodynamics of Landau and Lifshitz. 
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I. INTRODUCTION 

The relativistic generalization of the hydrodynamic 
theory, as developed by Landau and Lifshitz Q, leads 
to differential equations of the parabolic type. Consider 
for instance the diffusion equation for the viscous shear 
flow. In a frame in which the background equilibrium 
state is at rest ( "comoving frame" ) , it is 



[c- 2 {e+p)d t -r]d 2 x \ Su l (x, t) = , 



(1) 



where Su l with i = y, z is the perturbation in the 
transverse velocity. Equation ([!]) is the archetype of a 
parabolic equation. In non-relativistic physics, it pro- 
vides an excellent description of a wide range of phys- 
ical phenomena, as countless experiments have shown. 
Within the relativistic framework, however, it seems 
to fail, as two deficiencies become evident: The first, 
acausality, refers to the fact that Eq. (Q) allows for prop- 
agation of signals with arbitrarily large velocities. (This 
should be, and indeed has already been, a worry in the 
Galilean hydrodynamics, since one has very definite ideas 
of the velocities of the constituent microscopic particles, 
which represent an upper limit of the signal velocity in a 
dilute system of hard-core interaction.) The second de- 
fect is an instability found by Hiscock and Lindblom , 
who showed that Eq. (Q) develops a solution that grows 
exponentially with time in any non-comoving Lorentz 
frames. It is noteworthy that the growth time scale was 
found to be microscopically short. 

To overcome these deficiencies, extended fluid theories 
were put forward which start from the hydrodynamic the- 
ory but include additional dynamic variables; see e.g. ||. 
The resultant larger set of phenomenological coefficients 
can be chosen such that all the equations are hyperbolic, 
ensuring causality and stability. The price for this is 
twofold: A rather more complicated theory and the dif- 
ficulty of finding a universally accepted set of additional 



variables (except perhaps in dilute systems). In fact, re- 
cently it has been shown by Geroch || and Lindblom || 
that the complicated dynamical structure which ensures 
causality is unobservable. The evolution of any physi- 
cal fluid state according to any causal theory results in 
energy-momentum tensors and particle currents that are 
experimentally indistinguishable from the respective hy- 
drodynamic expressions. 

In this paper, we take a step back and again fo- 
cus on the simpler and more universal parabolic equa- 
tion (0). While acausality is an expected feature of 
parabolic differential equations, already discussed in the 
non-relativistic context 0-0]) instability is not. Being 
absent in a comoving frame, one is astonished at its ap- 
pearance in non-comoving frames. Our main purpose 
thus is to examine the origin and physical relevance of 
the instability. 

The paper is organized as follows: In Sec. [lj some gen- 
eral mathematical aspects of partial differential equations 
are reviewed. We especially recall the intuitive meaning 
of the characteristics. Section III is devoted to the one- 
dimensional diffusion equation in a comoving frame. We 
discuss the problem of causality and examine the two 
types of Cauchy problems that can be formulated with 
respect to a parabolic equation. In Sec. IV the discus- 
sion is generalized to non-comoving Lorentz frames. We 
especially scrutinize the solution that Hiscock and Lind- 
blom have identified as an instability, and show that it is 
a result of an ill-posed initial value problem. In Sec. |v| 
we briefly examine the general case, in which the spatial 
dimension of the diffusion equation is greater than one. 



II. PRELIMINARIES 

Let us first review some general aspects of partial dif- 
ferential equations that can be found in standard text- 
books on mathematical physics (eg. Q). Consider a lin- 
ear partial differential equation of second order for the 
unknown function t?(x, t). It can be written, most gen- 
erally, in the form 

[Adl + 2 B 8 X d t + C d 2 } d + F(x, t, tf, d x #, d t d) = 0, 

(2) 

where A, B and C are given functions of the two inde- 
pendent variables x and t. Depending on the value of 
D = B 2 — AC at a given point, Eq. (||) is referred to be 
of the elliptic {D < 0), parabolic (D = 0), or hyperbolic 



1 



(D > 0) type in this point. In the following we shall 
restrict ourselves to the cases in which A, B and C are 
constants. The type of Eq. (0) then remains unchanged 
throughout the entire region. 
The equation 



A (d x (p) 2 + 2B (d x (p)(d t <p) + C (d t <p) 2 = 



(3) 



is called the equation of characteristics of the partial dif- 
ferential equation (ph. Correspondingly, the family of 
curves, 



(p(x, t) — const , 



(4) 



with ip(x, t) satisfying Eq. (g), is called the family of 
characteristics. We collect the following facts: 

(i) An equation of the hyperbolic type has two dis- 
tinct families of real characteristics, an equation of the 
parabolic type has only one; an elliptic equation does not 
have real characteristics. 

(ii) The equation of characteristics is invariant with 
respect to arbitrary transformations of the independent 
variables, x — x(x, t), t = t(x, t). This implies that, if 
ip(x, t) is a solution of Eq. (|3|), and if f{x, t) transforms 
into <p(x, t), then (p{x, t) is a solution of the equation of 
characteristics accompanying the transformed differential 
equation. 

(hi) The outer real characteristics that pass through a 
given point (xq, to) bound the domain of influence £Iq of 
this point. If we consider the variable t as the time and 
think of the solution t9(a;, t) as a quantity that varies in 
x-space with time t, then this means that for t > i the 
solution in the region outside Qo is not influenced by the 
initial data given at (xq, to). As an example, consider 
the telegraph equation 



a 



1 



-d 2 

2 °t 



r) = 0, 



(5) 



with given constants a, v > 0. Since D ~ a 2 /v 2 > 0, it 
belongs to the hyperbolic type. Its equation of character- 
istics is (d x (p) 2 — {dtf) 2 jv 2 = 0, which gives two families 
of characteristics, 



(p(x, t) — x ± v t — const. 



(0) 



Taking some point (xq, 0), the respective domain of in- 
fluence is thus bounded by x ± v t — xq , implying that 
the effects of the initial data propagate at finite velocity 
v. 

For equation (Q), the most general Cauchy problem is 
formulated in the following way: Let S be some smooth 
curve given in the space of the variables x, t. With each 
point (x, t) £ S there is associated some direction n not 
tangent to S. The problem now consists of finding, in 
some neighborhood of S (either on one or both sides 
of the curve), a solution {}(x, t) which satisfies the pre- 
scribed Cauchy conditions 



6\s = Q (x, t) 



dn 



&i{x, t). 



(7) 



It is important to note that the domain in which the un- 
known solution has to be determined is not specified be- 
forehand. So generally, the initial manifold 5* lies within 
the domain of definition of the solution. 

A problem is said to be well-posed if it has the follow- 
ing properties: The solution (i) exists, (ii) is uniquely de- 
termined, and (iii) depends continuously on the assigned 
data. The last requirement is imposed in connection with 
the fact that the initial data of physical problems are de- 
termined experimentally and so small errors occur. It 
is thus necessary to be sure that the solution does not 
depend essentially on the measurement errors of these 
data. 



III. THE DIFFUSION EQUATION IN A 
COMOVING FRAME 



A well-known example of Eq. 
parabolic diffusion equation 

dtd-adl$ = 0. 



is given by the 



where a > 0. With A — —a and B = C = its equation 
of characteristics (||) takes the form (d x tp) 2 — 0. Hence, 
the characteristics are given by the one-parameter family 
of lines 



</j(x, t) = t = const. 



(9) 



It is obvious from this that the diffusion equation allows 
for the propagation of disturbances with infinite velocity. 
Indeed, the initial value of ■& at the point (xq, to) has 
influence on the solution $(x, t) in the whole semi- infinite 
region t > to- 

This fact, however, does not in practice cause any com- 
plications: The superluminal propagation speeds are as- 
sociated only with variations on microscopically small 
time and length scales (and with amplitudes of the or- 
der of thermodynamic fluctuations). On these scales a 
macroscopic description loses validity, and the diffusion 
equation and its solutions break down accordingly. The 
acausal consequences are therefore precluded by restrict- 
ing the solutions of Eq. (^) to the hydrodynamic range 
of validity, 



1 



1 



(10) 



with $0 denoting the constant part of and £, r the 
characteristic distance and time between collisions of par- 
ticles, or elementary excitations. For gases this has been 
first demonstrated by Weymann || . 

As an example, consider the solution of Eq. (||) that 

satisfies the initial condition $|t = o = °do + {A/d) e~ x / 2d . 
It is 



d{x , t) = tf + {A/^d? + 2at) e- x 12 {d +2at \ (11) 



2 



which represents a Gaussian distribution with width 
2 at. As can be seen from the first of Eq. (10), 



this solution is of physical significance only in the inter- 
val 



|x| <€Ld 2 /£+{2a/£,)t. 



(12) 



The maximum speed at which measurable information is 
transmitted is thus of the order of u max ~ which 
is far less than the speed of light for typical values of 
the kinetic coefficient a and the microscopic length £. 
(Taking a as the heat conductivity, one has for iron at 



room temperature: a 10 m /s and £ 10 



So 



^max ~ 10 3 m/s, which is of the order of the speed of 
sound.) 

The most common initial-value problem for Eq. (^) is 
to find a solution $(x, t) for which the Cauchy condi- 
tions (^) are prescribed at the initial manifold t = 0. 
However, since the line t = is a characteristic here, 
the Cauchy data cannot be prescribed independently, but 
must satisfy a compatibility condition Q. For instance, 
if the Cauchy conditions are given by J?|t=o = ©o(x) 
and (9 t #)| t=0 = 6i(x), we get from Eq. (JsJ) , for t = 0, 
0i (x) = ad 2 Q>o(x). The characteristic Cauchy problem 
for the diffusion equation is therefore posed in the fol- 
lowing way: In the region t > find a solution "d(x, t) 
satisfying the initial condition 



0|*=o = e (x) 



(13) 



As is generally known |9|, this problem is well-posed 
for arbitrary, smooth functions 0o(x) that have a well- 
defined Fourier transform. In the special case of a pure 
exponential, Oo = Ae lkx (fc € H), the solution takes the 
form 



i?(x, t) = Ae 



ikx+Ft 



T(k) = 



(14) 



which is unique and stable for t > 0. 

Quite another situation arises when the Cauchy data 
are prescribed on a non-characteristic curve, say eg. x + 
bt = 0, with some finite 6 £ R. The diffusion equation 
then gives two modes, one bound and the other divergent 
with |x+6t| — > oo. A well-known example is the so-called 
"sideways problem": In the half-space x > (— oo < t < 
oo) we seek a function $(x, t), which satisfies Eq. (||), 
and which attains on the non-characteristic line x = 
the Cauchy conditions 



0|. 



Ae 1 ' 



(9xi?)|x=o=5e i 



(15) 



with w € M. As can be easily seen, this problem is 
generally solved by 

tf(x, t) = (Sdi e {l+1)AlX + (5i? 2 e ('+ 1 ) A ^) e wt , (16) 



Ai )2 (w) =±Vw/2o. (17) 
which represents a superposition of two modes 



# 1>2 (x, t) = «n(w) e *(±«v^7to+«t) e ±aiVS7to ( ( 18 ) 

with the amplitudes 2 (w) being determined by the 
two Cauchy conditions (Eq). One has 



l)y/2a/u. 



(19) 



So, without imposing further restrictions, the general so- 
lution ( |l6| ) explodes exponentially as x increases. More- 
over, it shows discontinuous dependence on the initial 
data (l7J: As ui — > oo, Eq. ( |l6|) is bounded on the ini- 
tial line x = but grows like exp(xy/u>/2a) for any 
x > 0. Consequently, the non- characteristic Cauchy 
problem ( plf ) is not well-posed. 

Recalling the intuitive concept of the characteristics 
(cf. Sec. y), the reason for this becomes obvious. The 
initial data given at some point (0, to) affect the value 
of the solution exactly in those points which lie in the 
domain of influence SIq = {(x, t)\x eE,(> to}. Solving 
Eqs. (||), ( p"5| ) thus yields the solution for x > as well 
as for x < 0. Confer the two modes of Eq. (|l8|). They 
describe damped waves which carry the initial data given 
at x — to the left and to the right. Physical intuition 
suggests, however, that $i and $2 do exist only for x < 
and x > 0, respectively. A superposition in the form of 
Eq. ( |l6| ) does not make any physical sense. 

Now, the sideways problem will not be well-posed un- 
less the behavior at infinity is prescribed. In fact, the 
physically realistic assumption that & be bounded as 
x — > 00 (TJ leads to a solution that exists, is unique, and 
depends continuously on the initial data given at x = 0. 
The correct formulation of the sideways problem thus is 
the following: Find the bounded solution of Eq. (||) in 
the region x > (—00 < t < 00), satisfying 



0| a =o = Ae* 



(20) 



Note that the boundedness condition "replaces" the sec- 
ond Cauchy condition. 

Equation (15) is an example of a non-characteristic 
Cauchy problem with the initial data given on a line that 
is timclikc. For the following, it is important to consider 
also the case where the Cauchy data are given on a space- 
like line. As we shall see, this problem is not well-posed 
either. 

In the half-space {fx, f)|x S K, ct — fix > 0} we seek 
the solutions of Eq. (g) satisfying periodic initial data on 
the non-characteristic spacelike line ct — (3x — (f3 £ 
H, < [3 < 1). Here c is the speed of light. With 
7 = (1 - [i 2 )- 1 ' 2 , x^ = (x, ct), if v = diag(l, -1), and 
= ~-f((3, 1), e M = 7(1, (3) respectively denoting the 
timelike and spacelike unit vector normal and parallel to 
the initial line ct — fix = n^Xu — 0, the general Ansatz 



§{x, t) 



k £ R, T gC (21) 



takes the form 
tf(x, t) 



ifc 7 (2;-/3ct)+r7(t-/3a:/c) 



(22) 
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with $\ c t-px=Q ~ exp(ifc7 l x). Inserting Eq. ( |22| ) into 
the diffusion equation (||) yields 



/3 2 



la^jT 2 - (l + 2i 7 a^fc ) T - 7 afc 2 + i/3cfc = . (23) 



For k the two roots r 12 (fc) are complex. The real 
parts are given by 



r R1 ,K 2 = \c ± J l -c 2 + y ^c 4 + lc3 7 -3 aA: 2 ; (24 ) 

where C = (ja(3 2 /c 2 )^ 1 . From this one finds the in- 
equalities 



r R1 + r R2 = oo 



(25) 



r R1 r R2 = \c 2 - J --c 4 + -c^ak 2 < o , (26) 

8 V 64 4 

which imply that exactly one mode grows exponentially 
with time. Now, as k — > oo, the general solution diverges 
like exp( 7A / c 3 /2a(3[t — (3x/c]\fk) for any t > (3x/c, while 
the initial values are bounded on J - (3x/c = 0. Hence 
small changes in the initial data would cause considerable 
changes in the solution, indicating that the problem is not 
well-posed. 



IV. THE DIFFUSION EQUATION IN A 
GENERAL LORENTZ FRAME 

Next we turn our attention to the diffusion equation 
in an inertial frame K, in which the medium moves with 
constant velocity v in the negative x-direction. Employ- 
ing the Lorentz transformation rules d t = 7& — 7 i | 9j and 
d x = jd x — 7 (u/c 2 )9j, we get from Eq. (g) the boosted 
diffusion equation 

2 

7 {d- t -vd^V-aj 2 (dl - 2 Jrd* d- t + ^df) = 0. 

(27) 

It is straightforward to verify that this equation still be- 
longs to the parabolic type. In a covariant language, it 
is rewritten as 



u^da <Q-a A^dn d v <& = . 



(28) 



where A^ = if v + c^u^u" , u* = j(-v, c), <9 M = 
(d s , d ci ), and rf = diag(l, -1). According to Eq. (§), 
the equation of characteristics becomes 



7 2 (d x <f~^ d~ t ip 



= 0. 



(29) 



so that the family of characteristics comes out as the 
general solution of the ordinary differential equation 
7 v dx + 7 c 2 dt = 0, or equivalcntly, 



u M dx* = . 



One finds 



<p[x, t) = J t + •y — x — t — const . 



(30) 



(31) 



with t being the proper time measured in a comoving 
frame K. Recalling Eq. ([)]), we explicitly see that the 
characteristics are invariant under Lorentz boosts; cf. 
Sec.JI 

With Eq. (|3|), the relativistic equivalent to the char- 
acteristic Cauchy problem (|l3|) is the following: Find a 
solution to Eq. (27) which satisfies prescribed values at 
the initial manifold t + (v/c 2 )x = 0. Clearly, this prob- 
lem is well-posed only for t+(v/c 2 )x > 0. Its exponential 
solution can easily be found from Eq. (|l4|) by making use 
of the Lorentz transformation, t = ft + j(v/c 2 )x and 
x — jx + jvt. 

Now, the solutions of Eq. (]27]) that Hiscock and Lind- 
blom H examine satisfy periodic initial data ~ e lkx 
(k £ K) on the non-characteristic line t = 0. Taking 



•Q{x, t) 



jkx+rt 



(32) 



they find the dispersion relation 



ic^-rf 2 - ( 1 + 2i"fa^rk] f - 7afc 2 + ivk = , (33) 
c 4 V c 1 ) 

see Eq. (67) of Ref. §. For k ^ the two roots f 1)2 (ft) 
are complex. The real parts satisfy the conditions 



>0, 



rRir R2 



c 2 fc . 

r tl < o 



(34) 



(35) 



which imply that one of the two modes diverges as 
t — > oo. ^From this the authors of Ref. |Q conclude that 
the Landau-Lifshitz theory of relativistic hydrodynam- 
ics is unstable (in the sense that small spatially bounded 
departures from equilibrium will diverge with time) and 
hence must be abandoned as a physically realistic theory. 

We disagree with this conclusion. The reason is that 
the above non-characteristic Cauchy problem is not well- 
posed, in complete analogy to its non-relativistic equiva- 
lent, Eqs. (|), (||). [Note that Eq. (§f) equals Eq. (||) 
if /3 = v/c] Its solution does not depend continuously 
on the initial data and therefore does not meet one of 
the three mathematical requirements to be posed with 
regard to realistic physical problems (cf. Sec. |l|). So 
it is the type of the initial-value problem here which is 
physically unacceptable, and not the instability of the 
resulting solution. 

Nevertheless, it is illuminating to interpret the 
Hiscock-Lindblom solution in terms of wave propagation 
(while ignoring the fact that the domain of definition of 



4 



the above Cauchy problem is restricted to the half space 
{(x, t)\x e IR, t > 0}). The general solution consists of 
two damped waves traveling to the regions t > and 
t < 0, respectively. Referring again to the discussion of 
the sideways problem, Sec. [II, this is a consequence to be 
expected. Since the initial manifold crosses the charac- 
teristics (|3~l|), the information (given at t = 0) propagates 
(along the characteristic lines t+ [v/c 2 ]x = const) in both 
directions. One thus gets two modes, each transporting 
the effects of the initial data in the respective region. 

It is important to note that the appearance of a mode 
running backwards in time in non-comoving frames is 
directly related to the infinite signal speeds accompany- 
ing the parabolic equation. In fact, it is well-known ]Tl| ] 
that "superluminal acausality" in a Lorentz frame K 
causes "chronological acausality" in any other Lorentz 
frames K. For example, consider a process in a frame 
K whereby an event V\ = (x, t) causes another event 
"Pi = {x + Ax, t + At) at superluminal velocity U > c. 
Let the time difference be At > so that Ax = U At > 0. 
Then, in a frame K , moving with velocity v relative to 
K, we have At = 7 At -7 (v 2 /c) Ax = 7 At (1-vU/c 2 ). 
If c 2 /U < v < c, this yields At < 0. This means that in 
K the signal goes backward in time, or equivalently, that 
the response Vi precedes the stimulus V\. 

Now, as we have seen in Sec. 



[II, the superluminal 



acausality in K does not lead to any physical conse- 
quences; it is automatically precluded by a restriction 
to the hydrodynamic range of validity. Accordingly, we 
see from Eqs. (|34|) and (|3^) that the damping rate of the 
mode which propagates in the past is bounded below by 
r R+ > c 4 /jav 2 . With t denoting the characteristic mi- 
croscopic time in K (referred to as the collision time in 
dilute systems) , and f = 7T the respective time in K, we 
thus have 



» 1 



(36) 



(For instance, taking a as the heat conductivity of a 



piece of iron, one has a ~ 10 m /s and r 



10 



so that r R+ f pa 10 .) This shows that the mode run- 
ning backwards in time t decays on a time-scale that is 
much shorter than the microscopic one, and hence is far 
outside the hydrodynamic regime. 

Let us now come back to the non-characteristic Cauchy 
problem. Since the domain of definition of the solution 
is {(x, t)\x &JR,t> 0}, the mode reflecting the chrono- 
logical acausality is defined only for positive times t. As 
a consequence, the superposed solution does not depend 
continuously on the initial data, and the problem is ill- 
posed. With other words, prescribing Cauchy conditions 
on a non-characteristic hypersurface in a 2-dimcnsional 
flat spacetime does not lead to a realistic physical prob- 
lem. Such a problem might be obtained only if one pre- 
scribes the value of the solution at infinity (cf. the side- 
ways problem in Sec. Ill), or if the Cauchy data are pre- 
scribed on one of the characteristics. 



Now, what arc the physical conclusions after the math- 
ematical ones have been drawn? Given the fact that 
the standard Cauchy problem is ill-posed in the Landau- 
Lifshitz theory, that one cannot arbitrarily prescribe ini- 
tial data of physical quantities on a boosted time slice, 
are we to conclude that the theory is defunct - and to 
be replaced by one in which the Cauchy problem is well- 
posed? Is the capability to accommodate arbitrary initial 
data on boosted time slices a sine qua none for a healthy, 
physical theory? We believe the answer is no, for two dif- 
ferent reasons, although we concede that an affirmative 
answer may also be upheld. Our reasons are 

• All differential equations belonging the the hydro- 
dynamic theory and accounting for coarse-grained, 
irreversible physics - including the diffusion equa- 
tion - possess a preferred inertial system, in which 
the material is at rest. This is not true with re- 
spect to microscopic theories for which the Cauchy 
problem is known to be well-posed in any inertial 
frame, consider for instance the vacuum Maxwell 
equations. 

• More importantly, in our simple example consid- 
ered above, we actually know what the initial data 
on a boosted time slice t — correspond to in the 
rest frame, namely to an ill-posed, sideway problem 
arising from data on a non-characteristic spacelike 
line ct — "f(v/c)x = 0. This connection is com- 
pletely general, and especially independent of the 
Landau-Lifshitz theory. Therefore, an insistence on 
the well-posedness of Cauchy problems for boosted 
systems amounts to requiring a new type of differ- 
ential equation for the rest frame, ones for which 
the sideway problem is well-posed. Now, there 
are few equations in physics which are better con- 
firmed than the rest frame diffusion equations as we 
know them: for temperature, velocity and concen- 
tration. And it seems highly unlikely that one can 
change something as basic as the ill-posedness of 
the sideway problem without destroying the agree- 
ment with the experiments. 

Resigning to the fact that the standard Cauchy problem 
on a boosted time slice may be ill-posed, then clearly, 
there is no reason whatever to abandon the Landau- 
Lifshitz theory. 



V. GENERALIZATION TO CASES OF MORE 
THAN ONE SPATIAL DIMENSION 

Let us finally examine the general case in which the 
spatial dimension of the diffusion equation is greater than 
one. Consider the covariant equation (p8|), with u run- 



ning now from 1 to 4. As we have seen in Sec. IV, its 
characteristics are determined by the ordinary differen- 
tial equation 



5 



Ufj, dx^ = . 



(37) 



which assigns to each point x^ an infinitesimal character- 
istic 3-flat normal to the direction of the time-like velocity 
u^ at that point. The crucial point now is that generally 
these infinitesimal 3-flats do not integrate to 3-surfaces 
tpix 1 *) = const. 

To see this, let us suppose that a general integral 
^(x^) = const exists. Then, we have 



dip = {dn cp) dx>* = , 



(38) 



and hence 



dp ip — A Ufj, , (39) 
with some A = A(cc AI ). Since d^d v Lp = d v <p, this yields 
dfj. u v — d v ^ = A -1 (uft d v \ — u v <9 M A) , (40) 
from which e^ va ^u v d a up — or 
u a {duU^ - dfj,u v ) 

+U f _ l (8 a U v - 3 V U a ) + U v (dft U a - 8 a Uft) = (41) 

results. Equation ( |4l| ) is a necessary condition for intc- 
grability. It can be proved (see eg. 12 1) that it is also suf- 
ficient, ie., if it is satisfied, a general integral exists. [Note 
that Eq. ( fill ) is identically satisfied for a, /i, v = 1,2. 
Thus, in the case of one spatial dimension, characteris- 
tics can always be found.] 

Multiplying the condition of integrability with u a , one 
finds 



A^A^(daUf3-df3Ua) = 0- 



(42) 



So, characteristic hypersurfaces exist only if the back- 
ground equilibrium state is non-rotating. This fact has 
led various authors to believe that generally the parabolic 
diffusion equation lacks an initial- value formulation, and 
hence turns out not to be viable. This, however, would 
be true only if the non-characteristic Cauchy problem 
is ill-posed, just like it is in the case of one spatial di- 
mension. Unfortunately, very little seems to be known 
about the non-characteristic Cauchy problem in space- 
times with more than one spatial dimension. Further 
work should give clarity here. 
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